arXiv:1507.02775vl [math.AP] 10Jul2015 


A NEW FORMULA FOR THE ENERGY OF BULK SUPERCONDUCTIVITY 


AYMAN KACHMAR 


Abstract. The energy of a type II superconductor submitted to an external magnetic field of 
intensity close to the second critical field is given by the celebrated Abrikosov energy. If the 
external magnetic field is comparable to and below the second critical field, the energy is given 
by a reference function obtained as a special (thermodynamic) limit of a non-linear energy. In 
this note, we give a new formula for this reference energy. In particular, we obtain it as a special 
limit of a linear energy defined over configurations normalized in the L^-norm. 


1. Introduction 

1.1. A background. The behavior of a superconductor subjected to an external magnetic field 
varies as the intensity of the applied magnetic field changes. That has been observed early 
in the physics literature on theoretical and experimental grounds. There are two important 
key observations regarding a specific class of materials, called type II superconductors, namely, 
the formation of Abrikosov lattices and the persistence of surface superconductivity. Abrikosov 
lattices occur when the intensity of the magnetic held is near a special value, called the second 
critical field and denoted by Rc' 2 - If fh® intensity of the external magnetic held is increased 
above this value, then superconductivity disappears from the bulk of the sample and remains 
on a (part) of the surface of the material. This phenomenon persists until the intensity of the 
applied magnetic held reaches another special value, called the third critical field and denoted by 
Hcs ■ When the intensity of the applied magnetic held is increased further, superconductivity is 
destroyed everywhere in the sample, which is set into the normal state. The reader may consult 
the book of de Gennes [6] for the description of these important observations. Both phenomena, 
Abrikosov lattices and surface superconductivity where observed by theoretically investigating 
the Ginzburg-Landau model. Ginzburg and Landau proposed the model on a phenomenological 
basis to describe the response of a superconductor to an external magnetic held. 

Mathematically, the Ginzburg-Landau model is a functional dehned on a certain class of 
conhgurations. The physically relevant states of the superconductor are those corresponding 
to minimizing conhgurations (and critical points) of the functional. The Abrikosov lattice is 
distinguished by a special behavior of the minimizing conhgurations. The same applies for the 
surface superconductivity phenomenon. The two monographs mm contain many mathematical 
results regarding the Ginzburg-Landau model together with the discussion of their signihcance 
in physics. 

1.2. The Ginzburg-Landau model. Here we describe the Ginzburg-Landau model for a su¬ 
perconducting sample, occupying an inhnite cylindrical domain. The cross section of the cylinder 
is assumed a smooth and simply connected open subset Q of . 

The sample is subjected to an external magnetic held with direction parallel to the axis of the 
cylinder. The intensity of the external magnetic held is assumed constant. 

The superconducting material is distinguished by a characteristic parameter k > 0. When k 
is large, the material is of Type II. The intensity of the external magnetic held is denoted by a 
parameter hex- 

The behavior of the superconductor is described by a wave function ^ : H —)• C and a vec¬ 
tor field A : H —)• M^. The signihcance of and A is as follows. measures the density 
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of the superconducting Cooper pairs, whose presence is necessary to observe the superconduc¬ 
tivity phenomenon, and curl A measures the induced magnetic field in the sample, if present. 
At equilibrium, the configuration {'tp, A) minimizes the following energy, that we will call the 
Ginzburg-Landau energy: 

S{'ip,A) = j (^{V — iKHA)ip\^ +'^{1 — + K^H'^\cviT\A — 1\^ dx , (1.1) 

in the configurations space 

(V>, A) € C) X M^) , 

where denotes the usual Sobolev space. 

We introduce the ground state energy of the functional in dm as follows, 

Eg,(K,/7) = inf{f(V’,A) : A) eC) xM^)} , ( 1 . 2 ) 

Here we use the version of the functional as in [7] • In this version, the intensity of the external 
magnetic field is measured by /lex = nH, where > 0 is the the parameter marking the 
variation of the magnetic field. One way to detect the response of the superconductor to the 
external magnetic field is to mark changes in the energy as the parameter H changes. A large 
part of the mathematical literature is devoted to the computation of (11.21) when /lex is a function 
of K and K —)• 00 . The reader is referred to the two monographs mm for a detailed discussion on 
the behavior of the energy in (II.2p . In two special regimes, the energy in (II.2p is given to leading 
order by the Abrikosov and bulk constants introduced below in (II.6p and (EH respectively. 

1.3. The Abrikosov energy. We introduce the Abrikosov energy in a simple situation. More 
general situations are discussed in [1] but in the asymptotic limit considered here, they give 
rise to the same Abrikosov constant in (II.6p below. Let ii > 0 and consider the lattice in 
generated by the square 

Qr = {-R/2, R/2) X {-R/2, R/2 ). (1.3) 

Let us suppose that R? G 27rN. Denote by Pr the self-adjoint operator 

PR = -{y - iA^f in L^ag,per(Qi?) > 
defined via the closed quadratic form 

qR{u) = / \{V - iAo)uf dx . 

Qr 

Here, the vector potential Aq is defined as follows, 

Ao(xi,X 2) =-^(-X2,Xi) , (1.4) 

and generates a unit constant magnetic field, 

curl Aq = 1. 

The space b) and the form domain D{qR) are defined as follows, 

^Lg,per(< 5 R) = {u£ : ^(3^1 + R1X2) = e*^^i/^u(xi, 3 : 2 ) 

and n(xi, X2 + R) = X2)} , 

and 

D{qR) = ^^iag,per(QR) - G ■ (V - iAo)u G LL(]R2)} . 

The space peAQ r) a Hilbert space with the following inner product 

^Qr 

The spectrum of the operator is explicitly given by the Landau levels, 

(t(Pr) = {(2n - 1) : n G N} , 
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and, as long as B? G 27rN, all the eigenvalues have finite multiplicity. 

Now, we introduce the Abrikosov energy functional (in the square lattice). 


£Ah{u) 



(1.5) 


defined for configurations u in the finite dimensional space Er := Ker(Pfj — Id), the first 
eigenspace of the operator Pr. Minimizing the Abrikosov energy functional leads us to introduce 
the following quantity, 

eAb(-R) = inf{£:Ab(ii) : u G Er} . 

It is a known fact that [HlH]: 

R^G27rN 

Here, ^Ab is a universal constant that we will call the Abrikosov constant (or energy). It is 
known that ^Ab G [~ 51 0) ■ 

1.4. The reference ‘bulk’ energy. 

In this section, we recall how one can define the Abrikosov constant in (II. 6 p via a non-linear 
energy. Let b G (0,1] be a fixed constant, P > 0 and Qr be the square in (II.3p . Define the 
following non-linear functional in H^{Qr), 

£uk{u',b,R) = j ^6|(V - fAo)Mp - Irtp-h dx. (1.7) 

Here Aq is the magnetic potential in (11.411 . By minimizing this functional in various spaces, we 
get the following ground state energies, 

e^{b,R) = mf{£uk{u-,b,R) : u e H^{Qr)} , 

e^{b,R) =mi{£uk{u;b,R) : u £ H^{Qr)} , 

eP{b, R) = inf{£:bik(ii; b,R) : u£ Pmag,per(<3R)} • 

It is a known fact that [HEiiH] 

hm = Pblk(^) io£{D,N,p}), (1.8) 

R->-oo R^ 

where Pblk(') is a continuous and increasing function such that £'blk( 0 ) = —^ and Pblk(l) = 0 . 
The function £'blk(') is independent of the boundary condition and will be called the reference 
bulk energy. 

The Abrikosov constant in (|1.6p can be defined in the alternative way fflE], 


-^Ab 


-Eblk(&) 

lim 7---7 

(6 - 1)2 


(1.9) 


This formula displays a relationship between the non-linear simplified Ginzburg-Landau energy 
in (II.7p and the Abrikosov energy in (11.51) . 


1.5. The connection with the full GL functional. To illustrate how the quantities discussed 
so far are useful, let us cite the following two results from |15l IIP) . (We will use the following 
notation: If a(«:) and 6 (k) are two non-negative functions of k, then by writing a(K) <C 6 (k) we 
mean that a(«:) = 5{K)b{K) and 5{k) —)• 0 as «; —)• 00 .) 

(1) If 6 G (0,1] is a constant, H = bn and k -A 00 , then the ground state energy in (II.2p 
satisfies, 

H) = ^\n\Eukib) + o{k^) . 

(2) U H = K — h{k) and ^/K pl{k) <C k, then as k ^ 00 , 

Eg«(«:, H) = [k- P]2|L!|PAb + o{[k - Hf) . 


( 1 . 10 ) 
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1.6. The new formula. Here, we will define the function £'bik(-) via a non-linear eigenvalue 
problem. Let again 6 G (0,1], i? > 0, Qr be the square in (11.31) and Aq be the magnetic potential 
in dm). Let us define the linear functional, 


£iin{u;b,R) = / ^6|(V — iAo)u|^ — dx . 


( 1 . 11 ) 


' Qr 

We will minimize this functional in various spaces but for the constrained configurations 


/ 

•^Qr 


Ittr dx = 1. 


That way, we get the following ground state energies. 


m'^( 6 , R) = inf 


m 


m' 


N 


( 6 , R) = inf I 


\h,R) = inf I 


7 

/ \u\^ dx J 

'Qr / 

b, R^ 

ur dx 
Qr 

^lin(^) b, d?) 


: ueH^{QR)\{0}], 

: u €\ { 0 } } , 


lul'^ dx 


We will prove that, 

and that 
More precisely: 


lim 


'Qr 
m°(6, R) 


1/2 


: ueH, 


R—>-oo R 


— Enew (6) {oe{D,N,p}), 


Enewib) = --\/-2£'blk(^) • 


( 1 . 12 ) 


(1.13) 


7ag.per(Qii) \ {0} } . (l-M) 


(1.15) 


Theorem 1.1. Let b G (0,1). There exist two constants C > 0 and Rq > 1 such that, for all 
R> Rq and o G {D,N,p}, 

-(-2i?bik(&))'^'-^(-2i?bik(6))“'^' < < _(^_2Eu^{b))^/^ + ^. 

In particular, we may define the function £(blk(') formula, 

£b,t((.) = -i(/im (6 €(0,1)), 

with o G {D, N,p}. 

For the Dirichlet boundary condition, the lower bound is as follows 

and is valid for all b G (0,1] and R> 1. 


Let us compare the various ground state energies discussed so far. The Abrikosov functional 
in m is defined via a simple expression that does not involve differentiation operations but is 
minimized in the non-trivial space of the ground state eigenfunctions of the operator Pr. Among 
the remaining functionals we discussed, the expression of the non-linear functional in (|1.7p is the 
most complicated, but it is minimized in the space of Sobolev functions (this space is less com¬ 
plicated than the space of the ground state eigenfunctions of the operator Pr). The expression 
of the functional in ( 11 . 111 ) is linear, but again this functional is minimized for constrained config¬ 
urations. The three functionals serve in defining the Abrikosov constant L^Ab in (I1-6I) . but each 
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time one singles a simpler expression of the functional, a price is paid through a constraint in 
the definition of the ground state energy. 

Recently, there is a progress in the analysis of semi-classical non-linear eigenvalue problems 
with a magnetic field (cf. |12|1. The result in Theorem 1 1.1 1 mav fall in this area as well. 

The rest of the paper is decomposed into three sections. The proof of Theorem 11.11 occupies 
Sections [2] and [3J In Section HI we apply Theorem 11.11 to write a new proof of an important 
theorem by Ahnog devoted to the full Ginzburg-Landau functional (cf. [21 Thm. 3.3]). 


2. Proof of Theorem 11.11 Upper bound 


Let b G (0,1). This section is devoted to the proof of the following inequality 

m°(6,R) < -R(-2 ^bik(fo))'^" + C', (2.1) 

valid for and o G {D, N,p} and R > Rq, where G > 0 and Rq > 1 a-r® two constants that depend 
on b. 

Along the proof of f|2.1D . the following two lemmas are needed. 


Lemma 2.1. (]U] and |l]j 

There exists a constant C > 0 such that, for all 6 G (0,1], R> 1 and o G {D,N,p}, 

EuUb) - I < < ^bik(6) + ^ . 

For the Dirichlet boundary condition, the lower bound is 

e^{b,R) 


Ehik{b) < 


R2 


Lemma 2.2. There exists a constant G > 0 such that, for all b G (0,1], ii > 1 and o G {D, N,p}, 
if Ub^R is a minimizer of e°{b,R), then, 

—2i?^il/bik(b) — CR < / \ub^R\'^ dx < —2i?^il/bik(^) + Gii. 

J Qr 


Proof. The minimizer Ub^R satisfies the following equation 

- 6 (V — iAofub^R = (1 — |rt 6 ,_Rp)Rfe,i? in Qr , 

with adequate boundary conditions along the boundary of Qr (Dirichlet for o = D, Neumann 
for o = N and magnetic periodic for o = p). 

Multiplying the equation of Ub^R by Ub^R, integating over Qr then applying an integration by 
parts, we obtain after a rearrangement of the terms, 

I \ub,R\^ dx = £uk{ub,R) = e°{b,R). 

^ JQr 

Now, applying Lemma [2.II to estimate e°{b,R), we get the conclusion in Lemma [2.21 □ 


Proof of (|2.ip . Let 6 G (0,1], R > 1 and Ub^R be a minimizer of e°{b, R) for o G {D, N,p}. The 
ground state energy e°{b, R) is displayed right after introducing the bulk functional Tbik in (|1.7I) . 
We write using in particular the definition of m°{b,R), 


e°(b, R) — £h\k{ub,R) 

>m°{b,R)([ \ub,R\'^dx') +1/ \ub,R\'^ dx. 

\JQr / ^ JQr 


'Qr / ^ JQr 

Next, we estimate the L'^-norm of Ub^R by using Lemma [2.2l to write, for some universal constant 
G> 0, 

e°{b, R) > m°{b, R) ( - 2R^ Euk{b) - Gr)^^' + ^ - 2 R2 Euk{b) - CR] 
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Now, we estimate e°{b,R) as in Lemma l2.ll arrange the terms and get for a possibly new value 
of the constant C > 0, 

1 /2 

R^Eu\,{h) + CR>r(^ -2Euvih) - CR~^^ ^ m°{b,R) - i?2^bik(&) - CR. 


Now we select Rq sufficiently large (depending on b) such that the term —2£'blk(^) 
is always positive for R > Rq, then we divide both sides by R^ — 2Liblk(^) ~ CR~^ 
rearrange the terms above to get (EH). 


- CR-^ 
and 
□ 


3. Proof of Theorem O Lower bound 
This section is devoted to the proof of the following inequality 

m°{b,R) > -i?(-2^bik(6))'^'-C'(-2L;bik(6))”'^', (3.1) 

valid for some universal constant C > 0 and for all b € (0,1), R > 1 and o E {D, N,p}. 

Proof of (j3.ip . Let b E (0,1), i? > 1 and Wb^p> be a minimizer of m°{b,R) for o E {D,N,p}. 
The ground state energy m°(6, R) is displayed right after introducing the bulk functional £\i^ in 

dnn). 

Let us normalize in as follows, 




( - 2i?2^bik(6)) 


1/4 


-Wb,R ■ 


||m,R||L4(Ofl) 

The L^-norm of the normalized function satisfies 

WwIrWlRq^) = [-‘2R^Eu^{b)f/\ 
By definition of the functional in (|l.lip and m°{b,R), we see that. 


m°{b,R) = 


£rm{wb,R) 


= ( - 2R^Eu^{b)) ^>^£UwIr) • 


\Wb,R\\LA^Q^) 

Now, we write using in particular the normalization of ^ and the definition of e°{h, R), 


(3.2) 


= £uk{wlR) - 1 / \wlR\‘^dx 
^ JQr 

= ^blk(^fe,R) + E\^i]^{b) 

>e°{b,R) + R^EQi^{b). 


We estimate e°{b,R) from below using Lemma [2.11 to obtain, 

^lin«fl)>2i?Xlk(?^)-C’i?. 

We insert this into (|3.2p to get the inequality in (|3.ip . The improved lower bound for the Dirichlet 
boundary condition holds since (cf. Lemma l2.ll) 

e^{b,R)>R^Euk{b). 

Also, this last bound trivially holds for 6=1, since iiibik(l) = 0 and the spectral theory of the 
magnetic Laplacian with a Dirichlet boundary condition yields that e^{l,R) >0. □ 
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4. Application: Almog’s L'^-bound 

In [2], Almog estimates the L'^-norm of the Ginzburg-Landau order parameter for the three 
dimensional functional (the proof is valid for the two dimensional functional as well). This bound 
is of particular importance to estimate the error terms when seeking a fine approximation of the 
ground state energy, as in lain]- 

Here, we consider the same question as in [5] but in two dimensions. The method we give 
works in three dimensions as well, but we restrict to two dimensions for the sake of simplicity. 

We prove the following theorem: 


Theorem 4.1. Let A E (0,1). There exist two constants C > 0 and kq such that, if k > kq, 
Ak < H < k, and (i/i, A)k,j 7 is a critical point of m, then 



(4.1) 


where C Ll is any square of side-length 2 k and satisfying 

Qk C {dist(x,cin) > 


We stress again that the estimate in (14.ip is proved in [2] for three dimensional domains and 
when = Ll. The proof we give to Theorem 14.11 is based on Theorem 11.11 and differs from the 
one used in [2]- 

Sharper versions of the bound in (|4.1I) are given in [lOl 113) . However, these improved versions 
of (|4.ip are based on energy expansions of the form in (|l.inp . The proof of (ll.lflp requires a 
rough control of the order parameter similar to the one in (|4.1li . In |10) . a key element in the 
proof of (11.1011 was a strong L°° bound on the order parameter, namely (cf. [8l 110) 1 

C / H \ 

- ^172 + 

where w CC 12. Having ()4.ip in hand, we can derive the energy expansion in (ll.lOp without 
using the L°° bound in (14.21) . 

Unlike the hard proof of (14.21) . the proof we give to ()4.ip seems quite general, does not require 
too much regularity of the critical points and works for functionals having a similar structure as 
that in dni, e.g. the functional with a variable magnetic field or with a pinning term (cf. [31I5])- 

The rest of this section is devoted to the proof of Theorem 14.11 We will use C to denote 
positive constants independent of k and H. The value of C might change from one formula to 
another without explicit notice. 


4.1. Preliminaries. A critical point (fj, A) of the functional in (II.ip is a solution of the following 
system of PDE: 

r —{V — iKHA)‘^'ip = k‘^{1 — , 

< —V"*-curl A = (Ki/)~ilm('!/l (V — iKi/A)i/:), in 12 , (4-3) 

V ■ iy — iKH A)^} = Q , curlA = l on cA2 . 

We collect useful a priori estimates in: 

Lemma 4.2. Let A E (0,1). There exist two constants C > 0 and kq > 0 such that, if k > kq, 
Ak < H < k, and (i/i, A)K,i 7 is a solution of (14.31) . then, 

llV’lloo < 1 , 

and 

llcurlA- l||^i(j^) < Ck-^ . 

We refer the reader to for the proof of Lemma 14.21 Based on the estimates in Lemma 14.21 
one can construct the gauge transformation given in the next lemma (cf. |10l Eq. (5.30)]): 












8 


AYMAN KACHMAR 


Lemma 4.3. Let A G (0,1). There exist two constants C > 0 and kq > 0 such that the following 
is true. 

Suppose that k > kq, Ak < H < k, and ('0, is a solution of (14. 3p . Let I G (0,1) and 

Bi C Ll be a disk of radius 2£. There exists a function f G such that, 

V X G , I A(x) — (Ao(x) — V(/>(x)) I < CK~^i . (4.4) 

Here, Aq is the magnetic potential in (m. 

4.2. Local estimates - Proof of Theorem 14.11 Here we work under the assumptions in 
Theorem im We will estimate the following local energy of the critical configuration (fj, A): 


£< 


[ (\{V-iKHA)^|;f 

Qk V 




. .,4 


dx. 


(4.5) 


Let Q 2 k be the square having the same center as but with side-length i.e. twice 

the side-length of Q^- Obviously the square Q 2 k contains Let / G C^(Q 2 k) be a cut-off 
function satisfying, for all At > 1, 

/ = 1 in 0 < / < 1 and |V/| < in Q 2 k , 

where C is a constant independent of k. 

An integration by parts and the first equation in (|4.3I) yield (cf. |10l Eq. (6.2)]), 


£o,2K{f'ilJ,A) = Hi^f f (-l + ]-fA\'il;\'^dx+ [ iV/PiV’pdx 

d Q2 k \ ^ / J Q2k 


(4.6) 


<C. 


Note that we dropped the term involving —l + ^f‘^ since 0 < / < 1. The term involving |V/| is 
estimated using the bounds ||V’||oo < Ij |V/| < and \Q 2 k\ ^ Ck~^. 

Now we estimate the linear energy 


Ci 


,k(/V’,A)= [ (\{V - iKHA)fijj\‘^ - K^\fijj\‘^'] dx . 


(4.7) 


Choose C > 0 large enough such that Q 2 k, C B(j,^-i/ 2 , then apply Lemma 14.31 in B^,,-i /2 to 
get a function (j) G H‘^{Bq,^-i/ 2 ) satisfying the estimate in (14.4p . Notice that, using the Gauge 
invariance then the Cauchy-Schwarz inequality, 


^o,k(/V', A) = £o.« 




> 


[ ((1 - AC"^/^)|(V - mH'Ao)/V’e““^'^|^ - - CK^^^lfipl^^dx. 

JQ2. ^ ' 

Let 6 = (1 — R = kH = '/H and x^ be the center of the square Q 2 k- Apply 

the change of variable y = y/ kH {x — x^) to get, 

£o,.(/V’, A) > - CK^/^wmi, 

where m^{b,R) is the energy introduced in (|1.12p . We use the lower bound for m^{b,R) in 
Theorem o and Holder’s inequality for the term ||/a />||2 to get, 

£o,«(/V’, A) > ac3/2^-i/ 2 _ 2Ebik(fe))'/")i?||/V’lll - CK\\ml. 

Recall that R = '/H and insert the result into (|4.7I) and the right side of (14.61) to get, after a 
rearrangement of the terms, 

«-4/2 . 




+ 


l<c. 


(4.8) 


Two cases may occur: 
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Case I: 


Case II: 


2 

. 1/2 


(-|2Ebik(^)r/'-CK-V2) + 

Clearly, in both cases, (|4.8I) yields the following upper bound: 




5 




Since / = 1 in Qk, then (|4.9I) says that 


f f <2 k ^'^^|2£lbik(fe)|^^^ + 

\Jq^ ) 


(4.9) 


(4.10) 


Recall that h = (1 — i/^)— < 0. The assumption in Theorem 14.11 and the formula in (11.91) 

together yield that 

H 

- 1 


|2i?bik(?^)| <C|6-1| <C 


+ CK-V2. 


Inserting this into (|4.10p and remembering that \Qk\ = 2k we get the estimate in (|4.ip . This 
finishes the proof of Theorem 14.11 
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